Effective contact pairing forces from realistic calculations in infinite homogeneous 

nuclear matter 



O 






> 

m 

00 

p 
o 

O 



X 



N. Chamel 

Institut d'Astronomie et d'Astrophysique, CP-226, 

Universite Libre de Bruxelles, 1050 Brussels, Belgium 

(Dated: July 7, 2010) 

Non-empirical efTective contact pairing forces to be used in self-consistent mean-field calculations 
are presented. These pairing forces, constructed so as to reproduce exactly any given microscopic 
pairing gaps in infinite homogeneous nuclear matter for any isospin asymmetry, are given in ana- 
lytical form. As a by-product, this work provides an analytical solution of the BCS gap equations 
which could be applied to describe various many-body systems. 
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I. INTRODUCTION 

Self-consistent mean-field calculations using effective 
interactions have been very successful in describing the 
properties and the dynamics of a wide range of nuclei [1| . 
One of the most popular effective interactions are zero- 
range interactions of the Skyrme type because of the fast 
numerical computations which can thus be performed [^j ■ 
Even though Negele and Vautherin Q showed a long 
time ago how to obtain effective interactions using nu- 
clear many-body methods, a more empirical approach 
has been usually followed. A specific parametric form of 
the effective interaction is postulated and the unknown 
parameters are determined a posteriori so as to repro- 
duce a set of nuclear data selected according to a specific 
purpose. The non-uniqueness of the fitting procedure has 
thus lead to a large number of different parametrizations. 
Some of them may yield very different predictions when 
applied outside the domain where they were fitted [3l- 

This situation is particularly unsatisfactory for nuclear 
astrophysical applications which require the knowledge of 
nuclear masses for nuclei so neutron rich that there is no 
hope of measuring them in the foreseeable future; such 
nuclei play a vital role in the r-process of nucleosynthe- 
sis [5] and are also found in the outer crust of neutron 
stars J3J . Extrapolations far beyond the neutron drip line 
are required for the description of the inner crust of neu- 
tron stars where nuclear clusters are embedded in a sea 
of superfluid neutrons iT] . Such environments cannot be 
reproduced in the laboratory. However many astrophys- 
ical phenomena observed in neutron stars are precisely 
related to the physics of these outer layers @] . The need 
for more reliable extrapolations of these nuclear models 
has motivated recent efforts to construct non-empirical 
effective interactions and more generally microscopic nu- 
clear energy density functionals [9]. Unfortunately such 
ah initio nuclear energy density functionals able to repro- 
duce existing experimental nuclear data with the same 
degree of accuracy as phenomenological interactions are 
not yet available. 

For the time being, a reasonable approach is to fit 
these phenomenological interactions with as many avail- 



nuclear energy density functional 



able nuclear data as possible, not only from experiments 
but also from realistic calculations in infinite homoge- 
neous nuclear matter. Very accurate nuclear mass mod- 
els based on the Hartree-Fock-Bogoliubov (HFB) method 
with zero- and finite-range effective forces have thus been 
recently developed [lOl . Ill| . In particular in our model 
HFB- 17 [13, [l3, we have achieved our best fit ever to 
essentially all the available experimental mass data, the 
rms deviation for the set of 2149 measured masses of 
nuclei with A^ and Z > 8 lij being only 0.581 MeV. Sig- 
nificant improvements compared to our previous models 
have been made possible by a better treatment of pairing 
correlations first introduced in Ref. [i5[ and extended in 
Ref. [ij, [lj| . The effective interaction used in the pair- 
ing channel was constrained to reproduce the ^Sq pair- 
ing gaps as obtained by microscopic calculations using 
realistic nucleon-nucleon interactions. We have recently 
applied these effective forces to study the properties of 
the neutron superfluid phase in the inner crust of neutron 
stars [iq . 

The effective pairing interaction associated with our 
nuclear mass models requires the evaluation of one di- 
mensional integrals. Although its numerical implemen- 
tation is straightforward, using such interaction could be- 
come computationally costly for large scale calculations 
such as the determination of fission barriers. In this pa- 
per, a more tractable expression of the pairing interaction 
is presented by calculating the underlying integrals an- 
alytically. We will consider effective interactions of the 
Skyrme type in the particle-hole channel. But results can 
be easily generalized to other kinds of interactions. 



II. EFFECTIVE DENSITY-DEPENDENT 
CONTACT PAIRING FORCE 

The pairing interaction that we consider here acts only 
between nucleons of the same charge state q (q = n or p 
for neutron or proton respectively) and is given by 



pair 



in.rj) = -(1 - P,)v^^p,,ir),pp{r)] 5{nj) , (1) 



where Pa is the two-body spin-exchange operator, rij = 
Tj — Vj and r = {ri + rj)/2. Due to the zero range of 
this interaction, a cutoff has to be introduced in the gap 
equations in order to avoid divergences (for a review of 



the various prescriptions, see for instance Ref. |17|). In 
this work we include ah single-particle states whose en- 
ergy lies below Ag-|-eA, where Ag is the chemical potential 
and £a is an energy cutoff. Although microscopic calcula- 
tions in semi- infinite nuclear matter sug gest that pairing 



in nuclei is a surface phenomenon [18|, the density de- 
pendence of the pairing force is still poorly known. It is 
generally assumed that v^'^[pnjPp] depends only on the 
isoscalar density p = Pn + Pp and is often parametrized 
as [H 



'[Pn,Pp]^V^ 






(2) 



where po is the nuclear saturation density while VJ^, rjq 
and aq are adjustable parameters. The superscript A 
on VA is to remind that the pairing strength depends 
very strongly on the cutoff used in the gap equations (in 
principle changing the cutoff modifies also the other pa- 
rameters but the effects are usually found to be small). 
Effective interactions with 77^ = {rjg = 1) have been 
traditionally refered as volume (surface) pairing. The 
parameters in Eq. ^ are usually fitted directly to ex- 
perimental data. The standard prescription is to adjust 
the value of the pairing strength V^ to the average gap 
in ^^'^Sn [20] • However this does not allow an unambigu- 
ous determination of the remaining parameters rjq and 
aq . Systematic studies of nuclei seem to favor a so-called 
mixed pairing with r]q ^ 0.5 and 1/2 < a, ^ 1 [2lll22|. 

Even though such forces have been widely applied in 
nuclear structure calculations with some success, they 
lack a direct connection with realistic nucleon-nucleon 
forces. The reliability of these forces far beyond the 
domain where they were fitted is therefore not guaran- 
teed. In particular pairing forces fitted only to finite 
nuclei generally yield unrealistic pairing gaps in infinite 
nuclear matter |15l . |23| , thus rendering their ap plic ation 
to neutron-star crusts unreliable. Garrido et al. 24] pro- 
posed to determine the parameters of the pairing interac- 
tion in Eq. ([2|) by fitting the ^Sq pairing gaps in infinite 
symmetric nuclear matter as obtained by the realistic 
Paris potential in the BCS approximation. The para- 
metric form Eq. © has been recently extended by Mar- 
gueron et al. [23] who introduced an isospin dependence 
in the pairing strength in order to reproduce the ^ 5*0 pair- 
ing gaps in both symmetric nuclear matter and pure neu- 
tron matter as obtained from Brueckner calculations l26| . 
This approach assumes that the pairing interaction be- 
tween two nuclcons inside a nucleus is locally the same 
as the pairing interaction between two nucleons in infi- 
nite uniform matter. This should not be confused with 
the so-called "local density approximation" which sup- 
poses that the pairing field A(r) appearing in the HFB 
equations (see for instance Appendix A of Ref. [15|) is 
locally the same in finite nuclei and in infinite nuclear 



matter. Even though the coupling to surface vibrations 
is expected to contribute to pairing |27t. a local pair- 
ing theory seems a reasonable first step |28l . ^^ . Finite 
size effects can be subsequently taken into account by 
introducing density gradient terms in the expression of 
the effective pairing force as suggested in Ref. [30] in the 
general framework of the density functional theory. Ef- 
fective pairing forces with density gradients have been 
implemented in Ref. |3l| . 

The drawback of using a phenomenological pairing 
force in order to fit microscopic pairing gaps in infinite 
uniform matter is that the functional form of the associ- 
ated pairing strength is not a priori known. In particu- 
lar different pairing gaps require different expressions of 
v"^ "^ [pn, Pp] which can hardly be guessed due to the highly 
non-linear character of pairing correlations [321 1331 ]. It 
should also be remarked that the parameters of the pair- 
ing force depend on the effective interaction used in the 
particle-hole channel since effective masses appear in the 
gap equations. Fitting the parameters of the effective in- 
teraction to essentially all experimental nuclear masses 
both in the particle-hole and the particle-particle chan- 
nels while simultaneously reproducing microscopic pair- 
ing gaps in infinite uniform matter would thus be an ex- 
tremely onerous numerical task. 



III. MICROSCOPICALLY DEDUCED 
CONTACT PAIRING FORCE 

Actually as shown in Refs. [ij, [ij, [l^, the density- 
and isospin dependence of the pairing strength can be 
calculated exactly for any given ^ Sq pairing-gap function 
Aq[pn,Pp\ by solving directly the HFB equations in infi- 
nite uniform matter yielding 



v'"'[Pn,Pp] 



8n' 



Iq{Pn,Pp) \2M*{pn,Pp) 



3/2 



where Mq{pn, Pp) is the nucleon effective mass and 



IqiPn,Pp) = I dS. 



V^ 



V(^-Mg)' + A,(p„,pp)2 



(3) 



(4) 



The expression of M*{pn, Pp) for Skyrme forces can be 
found for instance in Appendix A of Ref. [Ta] ■ In Eq. Q , 
Pq — Xq — Uq is a rcduccd chemical potential {Uq being the 
mean field potential) which can be obtained by imposing 
the conservation of the nucleon particle number density 




However it is usually a very good approximation to re- 
place pq in Eq. (|4]) by the Fermi energy 



Aq) ^ "• ^Fq 
'^ 2M* 



(6) 



where fcpg = (STr^pg)^''^ is the Fermi wave number. The 

approximation fiq ~ £|? holds provided A^ <^ e^ which 
is typically the case |34| . 

Using Eq. ^ instead of a phenomenological expression 
such as Eq ([2]) guarantees that in infinite uniform matter 
the microscopic pairing gaps Aq{pn,Pp) will be exactly 
reproduced. It also ensures that the pairing strength 
is suitably renormalized for any change in the energy 
cutoff. The price to be paid is the evaluation of the 
one dimensional integral, Eq (|3]), for the neutron den- 
sity Pnir) and proton density pp{r) at all points inside 
the nucleus. The computational cost becomes significant 
in fully self-consistent 2D and 3D calculations {&]. How- 
ever we will now show that these numerical integrations 
can be avoided by making use of the weak-coupling ap- 
proximation. 

The integrand in Eq. ([4]) is all the more peaked around 



-(■Z) 



(«) 



£_ — fiq ^ £p as Ag is small compared to £p . This 
suggests to expand the integrand in powers of Ag/sp 
leading to 



I +00 

-"S ^ V F Z^ "? 
n=0 



(7) 



~ (") 
where the dimensionless coefficients /g are defined by 



~ («) ^ (2n)! 

' (l-2n)(n!)2 




ea/A, 



Ma; 



r<?VA, VT+X^ 



(8) 

In the weak-coupling approximation Ag <C e^' and 
Ag <C ea, only the first coefficient in Eq. ([7]) is usually 



(9) 



retained, i.e. /„ 



This approximation is 



equivalent in taking a constant density of single-particle 

states in the gap equations. Calculting the integral Iq 

and keeping lowest order terms in Aq/sp and Aq/sA, 
Eq. dU thus yields 



-(0) 



log 






(9) 



Inserting Eq. © in Eq. ^ and solving for Ag leads to 
the familiar expression (note that u^ ' < 0) 



A(o) 



(9) 

Ep EABXp 



V'"lM*kFq 



(10) 



Even though the weak-coupling approximation provides 
good results in the case of conventional BCS supercon- 
ductivity (33, it is less accurate in the nuclear context 
because of the large number of states involved in the 
pairing mechanism. Nevertheless the higher-order coef- 
ficients in Eq. (J8]) can be easily evaluated. Calculating 

the integrals Iq , keeping as before lowest-order terms 
in Aq/sp and Aq/sA, and summing all coefficients, we 
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FIG. 1: Neutron pairing strength v^"[pn, pp], as defined by 
Eq.(|3]), vs nucleon number density p in symmetric nuclear 
matter using the effective interaction underlying the nuclear 
mass model HFB-17 |13.[l3l]. The symbols are the results of 
numerically integrating Eq. Q. The dashed and solid lines 
were obtained using either Eq. ((9} or Eq. p2|l respectively. 
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+00 



^/g^"^^2yTT^-2iog i(i + yTT^) 



n=l 



-4, (11) 



where y = ea/ef^ Adding Eqs. ^ and ([TT]) in Eq. 



leads to 



HtH^, 



(12) 



with 



A(a;) = log(16a;)+2Vl + x-2 log (l + \/l + x) -4 . (13) 

The highly-non linear character of the pairing phe- 
nomenon is evident in Eqs. ([3]) and ([T2|) . These equa- 
tions also show that the density dependence of the pairing 
strength is intimately related to the choice of the pairing 
cutoff. The contribution of the latter is entirely contained 
in the function A(x). It can also be seen from Eqs. ([3]) 
and (J12p that the pairing strength vanishes whenever the 
pairing gap Ag goes to zero at finite density, as expected. 
Note that as a by-product we have also obtained a 
more accurate expression of the pairing gap. Substituting 
Eq. dni) in Eq. © yields 



A(o)exp(±A(y))j/-i/^ 



(14) 



where y ~ ea/ej? . 

We have tested the validity of the analytical expression 
given by Eq. (J12p as compared to the numerical integra- 
tion of Eq. Q using the effective interactions underlying 
our IIFB-17 nuclear mass model [l^, ll^l- The pairing 
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FIG. 2: Same as Figure [T] but for pure neutron matter. 



force used in this model was adjusted so as to reproduce 
the ^So pairing gaps in both symmetric nuclear matter 
and pure neutron matter as obtained from Brueckner cal- 
culations using the Argonne T^18 nucleon-nucleon poten- 
tial [26|]. The pairing cutoff was fixed at ea — 16 MeV. 
As can be seen in Figures [T] and [U Eqs. ^ and p^ pro- 
vide a very good estimate of the pairing strength for all 
densities and different isospin asymmetries. The figures 
also show that keeping only the first term in Eq. ([7]) with 

Iq given by Eq. ([9]) , leads to fairly good results except 
at low densities. 

At very low densities p — > 0, both the pairing gap Ag 
and the chemical potential vanish so that the integral 
Iq is then simply given by 2^JFJi. It is easily seen from 
Eq. ([3]) that the pairing strength remains finite in this 
limit and is given by 



y'^«[p->0] = - 



47r2 / n" 



'^\2Mq 



3/2 



(15) 



assuming that M*{p ^- 0) = Mq. The weak-coupling 
expression, given by Eqs. ([3|) and p^ . tends to the ex- 
act limit as shown in Figure [31 whereas Eq. ^ leads to 
a divergence (the corresponding pairing strength lies far 
below the range shown in Figure [3]). The reason for this 
discrepancy lies in the underlying assumption of a con- 
stant density of s.p. states, which is strongly violated at 
low energies e where the density of states decreases as 

The value of the pairing cutoff ea can be a priori ar- 
bitrarily chosen. However it has been argued [2J, |3g 
that in the limit p — > 0, the pairing strength should 
coincide with the bare force in the ^So channel, which 
in turn is determined by the experimental ^So nucleon- 
nucleon phase shifts. According to Eq. (fT5|) . specifying 
v'"^[p — ?> 0] is equivalent to specifying the pairing cutoff 
Ea- As shown in Fig. 1 of Ref. [23, the optimal value of 
the cutoff is £a ~ 7 — 8 MeV (note that ea is half the 
cutoff used in Ref. [3a)). Using the experimental phase 
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FIG. 3: Neutron pairing strength ii'^"[p„, pp], as defined by 
Eq.Q, vs nucleon number density p in symmetric nuclear 
matter (SNM) and pure neutron matter (NeuM) using the 
effective interaction underlying the nuclear mass model HFB- 
17 [l3, [l3l- The symbols are the results of numerically inte- 
grating Eq. Q while the solid line was obtained from Eq. (|12|l . 



shifts would thus remove the only free remaining param- 
eter. On the other hand, choosing such a low cutoff can 
deteriorate the precision of the weak-coupling approxi- 
mation, since Eq. ([12]) was obtained assuming Ag ^ ea- 
We have therefore calculated the relative error between 
the pairing strength calculated numerically and that ob- 
tained from Eqs. ([3|) and ([T^ . for different values of the 
pairing cutoff. Results are shown in Figs. [4| and [5] for 
symmetric nuclear and pure neutron matter respectively. 
We have checked that the differences come solely from 
the weak-coupling approximation, and not from the nu- 
merical method used to solve Eq. Q. It can be seen 
that for the effective interactions underlying our HFB-17 
nuclear mass model [ij, [iJl , the weak-coupling approx- 
imation is quite accurate since the largest relative error 
is of a few percent only for ea = 8 MeV. The precision is 
higher in neutron matter than in symmetric matter be- 
cause the ^Sq pairing gap of neutron matter is smaller 
than that of symmetric matter. Note that the pairing 
gaps we adopted in our HFB-17 mass model are rather 
large (the maximum of the gap function Aq(p„,pp) is 
about 4.8 MeV in symmetric matter and 2.6 MeV in neu- 
tron matter). These gaps were obtained from Brueckner 
calculations including medium polarization effects but 
without self-energy corrections [2g|. When both effects 
are taken into account, the maximum pairing gaps are 
much lower [26] . The precision of the weak-coupling ap- 
proximation would thus be even better with such gaps. In 
any case, if the constraint of reproducing the ^Sq nucleon- 
nucleon phase shifts is released, one is free to adjust the 
cutoff so as to achieve a better precision, as illustrated in 
Figs [land [5] 

In the more general context of the nuclear energy den- 
sity functional theory, the HFB equations can be ob- 
tained by minimizing the energy with respect to the 
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where Pqir) is the so-called local pairing density. The 
pairing field appearing in the HFB equations is then de- 
fined by 



A,(r) 



9gpair(?') 



(17) 



Because fpair('') depends on the nucleon densities, it con- 
tributes also to the mean-field potentials (see for instance 
Appendix A of Ref. [ia|)- Local pairing functionals can 
be constructed from infinite nuclear-matter calculations 
using Eqs. dSj, (fT2|) and (ITS]) . The construction of non- 
empirical (and non-local) pairing functionals using real- 
istic nucleon-nucleon interactions have been discussed in 
Ref. [33] (see also Ref. l9l| for a review) . 



FIG. 4; Relative error between the neutron pairing strength 
'c'^"[pn, Pv\ calculated numerically and that obtained from the 
weak-coupling approximation, vs nucleon number density p. 
The different curves correspond to different values of the pair- 
ing cutoff £a indicated in the plot. Calculations were carried 
out in symmetric nuclear matter using the effective interac- 
tion underlying the nuclear mass model HFB-17 [12.[l^. 
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FIG. 5: Same as Fig. |4] but in neutron matter. 



IV. CONCLUSIONS 



In conclusion, we have shown how to construct non- 
empirical effective contact pairing forces (or equivalently 
local nuclear pairing energy density functionals) using 
any given ^Sq pairing-gap functions Aq{pn,Pp) obtained 
from microscopic calculations in infinite uniform nu- 
clear matter with realistic nucleon-nucleon potentials. 
The resulting parameter-free pairing forces, embodied 
in Eqs. ^, (IT^ and (IT^ . can be easily implemented 
in nuclear-structure calculations. These forces could be 
helpful for understanding the origin of pairing in finite 
nuclei, and could be used to estimate the importance of 
bulk contribution as compared to finite-size effects. Al- 
ternatively, the analytical solution of the BCS gap equa- 
tions which we have obtained in Eq. (1141) , can be applied 
to calculate the pairing gaps in nuclear matter for any 
given contact pairing force (provided the conditions for 
the weak-coupling approximation remains valid). Even 
though we have been interested in nuclear pairing, the 
present results could be easily adapted to other contexts. 



normal and pairing density matrices for a fixed average 
number of neutrons and protons [33] (the link between 
different formulations_of the HFB equations is discussed 
for instance in Ref. [l5|). All that is needed is therefore 
the specification of the nuclear energy density functional. 
The pairing energy density associated with an effective 
pairing force defined by Eq. ([1}, is given by 



fpair(r) = i Y. v^''[Pn{r),Pp(r)]p,irf 



(16) 



q^n^p 
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